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Abstract 

We present a new algorithm for computing motorcycle graphs that runs in 0(n 4//3+e ) 
time for any e > 0, improving on all previously known algorithms. The main application of 
this result is to computing the straight skeleton of a polygon. It allows us to compute the 
straight skeleton of a non-degenerate polygon with h holes in 0(n\/h + 1 log 2 n + n 4 ^ 3+s ) 
expected time. If all input coordinates are 0(logn)-bit rational numbers, we can compute 
the straight skeleton of a (possibly degenerate) polygon with h holes in 0{ri\Jh + 1 log 3 n) 
expected time. 

In particular, it means that we can compute the straight skeleton of a simple polygon 
in 0(n\og 3 n) expected time if all input coordinates are 0(logn)-bit rationals, while all 
previously known algorithms have worst-case running time ui(n 3 ' 2 ). 

1 Introduction 

The straight skeleton of a polygon P is a straight line graph embedded in P, formed by the traces 
of the vertices of P when it is shrunk, each edge moving at the same speed and remaining parallel 
to its original position. (See Figure [!}) It has been known since at least the 19th century; for 
instance, figures representing the straight skeleton can be found in the book by von Peschka [32J. 
Aichholzer et al. |5j gave the first efficient algorithms for computing the straight skeleton, 
and presented it as an alternative to the medial axis having only straight-line edges. The 
straight skeleton has found numerous applications in computer science, for instance to city model 
reconstruction [27] , architectural modeling [26] , polyhedral surface reconstruction [7J HU [29] , 
biomedical image processing [16] . It also has a direct application to CAD, as it allows to 
compute offset polygons [18J. The straight skeleton has become a standard tool in geometric 
computing, and thus fast and robust software has been developed to compute it [TO ] [25 | [30] . 
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(a) Input polygon (b) Straight skeleton (c) Offset polygon 



Figure 1: The straight skeleton of a polygon, and a corresponding offset polygon. 



(a) A motorcycle graph (b) Induced motorcycle graph 

Figure 2: The motorcycle graph of a set of four motorcycles (a), and the motorcycle graph 
induced by the polygon from Figure [TJ 

The complexity of straight skeleton computation, however, is still very much open. The 
previously best known algorithms were the 0(n 17 / 11+e )-time algorithm by Eppstein and Er- 
ickson |18| . and the 0(ra 3 / 2 log 2 ra)-time randomized algorithm by Cheng and Vigneron |15j . 
The only known lower bound is Q(nlogrt), by a reduction from sorting [20]. In this paper, we 
give new subquadratic algorithms for computing straight skeletons. In particular, if all input 
coordinates are 0(logra)-bit rational numbers, we give an 0(ny/h + 1 log 3 n)-time randomized 
algorithm for computing the straight skeleton of a polygon with h holes. It is the first near- linear 
time algorithm for computing the straight skeleton of a simple polygon. 

Eppstein and Erickson |18j introduced motorcycle graphs so as to model the main difficulty 
of straight skeleton computation. We are given a set of n motorcycles, each motorcycle having 
a starting point and a velocity. Each motorcycle moves at constant velocity until it reaches 
the track left by another motorcycle, in which case it crashes. The resulting graph is called 
a motorcycle graph. (See Figure [2^.) The motorcycle graph is a special case of the straight 
skeleton, where each motorcycle is modeled by a small and thin triangle. Conversely, a polygon 
induces a motorcycle graph, where each motorcycle starts at a reflex vertex and moves with 
the same velocity as this vertex moves during the shrinking process. (See Figure [2)3.) Cheng 
and Vigneron [15] showed that computing the straight skeleton of a non-degenerate polygon 
reduces to computing this induced motorcycle graph, and a lower envelope computation; Huber 
and Held extended this proof to degenerate cases [25]. The lower envelope computation can be 
done in 0{n\Jh + 1 log 2 n) expected time if P has h holes. 

Previously, the bottleneck of straight skeleton computation was the induced motorcycle 
graph computation. This is our main motivation for designing a faster motorcycle graph al- 
gorithm. In this paper, we give an algorithm for computing a motorcycle graph that runs in 
0(n 4 / 3+e ) time, for any e > 0, improving on all previously known algorithms. Here is a brief 
description of our algorithm. For each motorcycle, we maintain a tentative track, which may be 
longer than its actual track in the motorcycle graph. We also maintain a set of target points, 
which contains the endpoints of the tentative tracks that have been created earlier for this 
motorcycle, and that it has not reached yet. Initially, the tentative tracks are empty, and then 
we try to extend them one by one, all the way to the destination point. If two tentative tracks 
cross, we retract them, by roughly halving the number of possible crossing points on each of 
them. After performing this halving, the tentative tracks do not intersect, and we can safely 
move the motorcycle that reaches the end of its tentative track first. Then we try to extend the 
tentative track of this motorcycle to its next target point, and repeat the process. An example 
is given in Appendix |B| 

Apart from obtaining better time bounds for straight skeleton computation, there are at least 
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two other reasons for studying motorcycle graphs. First, Huber and Held |25| used the idea of 
computing the straight skeleton from its induced motorcycle graph to design and implement 
a practical straight skeleton algorithm. So it is important, even in practice, to get a better 
understanding of motorcycle graph computation. Another motivation for studying motorcycle 
graphs is a direct application to computer graphics, for quad mesh partitioning [19]. 

Some of our results make no particular assumptions on the input, but we also present a 
few results where we assume that the input coordinates are 0(logra)-bit rational numbers. We 
believe that this assumption is sufficient for most applications. For instance, in the applications 
mentioned above, it is hard to imagine that the input polygons would have features smaller 
than lnm, and size larger than 1000km, so 64-bit integers should be more than sufficient. 



1.1 Summary of our results and comparison with previous work 



The main novelty in this paper is our algorithm for computing motorcycle graphs (Section 2.1 ). 
This algorithm is essentially different from the two previous algorithms \15\ ITS] that both sim- 
ulate the construction in chronological order. Our algorithm, on the other hand, does not 
construct the motorcycle tracks in chronological order: It may move some motorcycle to its 
position at time t, and then later during the execution of the algorithm, move another motor- 
cycle to its position at an earlier time t' < t. (We give one such example in Appendix [b]) This 
answers an open question by Eppstein and Erickson [18} end of Section 5], who asked whether 
the running time can be improved by relaxing the chronology of the events. 

Our algorithm uses two auxiliary data structures, one for ray shooting, and another for 
halving queries. Given a query segment on the supporting line of a motorcycle, a halving 
query returns a splitting point on this segment such that there are roughly the same number of 
intersections with other supporting lines on both sides. (See Section [l.2| ) The implementation 
of these data structures in different settings lead to different time bounds. 

For all our results, we use the standard real-RAM model [31], that allows to perform arith- 
metic operations exactly on arbitrary real numbers. But for some of our results, we make the 
assumption that all input coordinates are 0(logn)-bit rational numbers. It has two advantages: 
It yields better time bounds, and allows us to handle the straight skeleton of degenerate poly- 
gons. This improvement comes from the fact that, for bounded precision input, two distinct 
crossing points between the supporting lines of two pairs of motorcycles are at distance 2~°( logn ) 
from each other. It allows us to use a simpler halving scheme: Instead of halving a segment 
according to the number of intersection points, we use the midpoint according to the Euclidean 



distance. (See Section 3.3 and 4.1 



Arbitrary precision input. For our first set of results, the input coordinates are arbitrary 
real numbers, on which we can perform exact arithmetic operations. In this case, our new algo- 
rithm computes a motorcycle graph in 0(n 4//3+e ) time (Theorem [5]). This improves on the two 
subquadratic algorithms that were known before: the 0(n 17 / 11+e )-time algorithm by Eppstein 
and Erickson |18] , which was first published in 1998, and the 0(n 3 / 2 log n)-time algorithm by 
Cheng and Vigneron [15], which first appeared in 2002. 

We also give, in Section 3.2, an 0{Cn log 2 (n) min(C, log n))-time algorithm for the case of 



C-oriented motorcycles, where the velocities take only C different directions. This improves on 
the algorithm by Eppstein and Erickson [18] . which runs in 0(n 4 / 3+£ ) time when C = 0(1). 

Our last result with arbitrary precision input is an 0(n^ 3+£ +n\/h + 1 log 2 n) expected time 
algorithm for computing the straight skeleton of a polygon with n vertices and h holes. (This 
result does not hold for a degenerate polygon where two reflex vertices may collide during the 
shrinking process, as in Figure pi) It improves on the algorithm by Cheng and Vigneron [15] 
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Figure 3: A degenerate polygon and its straight skeleton (left). Two reflex vertices collide during 
the shrinking process, and a new reflex vertex appears (middle). The induced motorcycle graph, 
where a new motorcycle appears when the two other crash (right). 



which runs in 0(n 3 / 2 \og(n)+n\/h + llog 2 n) expected time. It also improves on the 0(n 17 / 11+e ) 
time bound of the algorithm by Eppstein and Erickson [18] , but their algorithm is deterministic 
and applies to degenerate cases. 



Bounded precision input. The following results hold when all input coordinates are 0(log n)- 
bit rational numbers. There has been recent interest in studying computational geometry prob- 
lems under a bounded precision model (the word RAM), for instance the computation of Delau- 
nay triangulations, convex hulls, polygon triangulation and line segment intersections [5| 112] . 

that a motorcycle graph can be computed in 0(n log 3 n) time if 



We first show in Section 



3.3 



the motorcycles move within a simple polygon, starting from its boundary. The only other non- 
trivial cases where we know how to compute a motorcycle graph in near-linear time seem to be 
the case where all velocities have positive ^-coordinate, which can be solved in 0{n log re) time 
by plane sweep, the case of a constant number of different velocity vectors [18], or a constant 
number of directions (Section 3.2). 



Then in Section |4.2[ we show that the straight skeleton of a polygon with n vertices and h 
holes can be computed in 0(ny/h + 1 log 3 n) expected time. This result still holds in degenerate 
cases. So with bounded-precision input, and if we allow randomization, it improves on the 
0(n 17 / 11+e )-time algorithm by Eppstein and Erickson [18j . When h = o(n/ log 2 n), it also 
improves on the 0(n 3 / 2 log 2 re)-time algorithm by Cheng and Vigneron |15| . which cannot handle 
all degenerate cases. 

In particular, our algorithm runs in expected 0(n log 3 n) time when h = 0, so it is the 
first near-linear time algorithm for computing the straight skeleton of a simple polygon. The 
previously best known algorithms run in w(re 3 / 2 ) time in the worst case |15| 118]. 



1.2 Notation and preliminaries 

For any two points p, q, we denote by pq the line segment between p and q. Unless specified 
otherwise, pq is a closed segment. The relative interior of pq is the open segment pq \ {p, q}. 
We say that two segments cross if their relative interiors intersect. 

The motorcycles are numbered from 1 to n. Each motorcycle i has a starting point Sj, 
moves with constant velocity and has a destination point di that lies in the ray (sj,Uj). (See 
Figure [4^,.) When p € Sidi, we denote by r(i,p) the time when motorcycle i reaches p, so 
p = Si +T(i,p)vl. The supporting line i% of motorcycle i is the line through Sj with direction Vi. 
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(a) Input (b) Motorcycle graph 

Figure 4: The input to the motorcycle graph problem (a), and the resulting motorcycle graph 
(b). 

Each motorcycle i starts at Si at time 0, and moves at velocity v*i until it meets the track 
left by another motorcycle and crashes, or it reaches di and stops. So motorcycle i crashes if 
it reaches a point p such that r(i,p) i~(j,p), for some motorcycle j that has not crashed or 
stopped earlier than r(j,p). If motorcycle i crashes, we denote by Ki the point where it crashes, 
called the crashing point. (See Figure ) Otherwise, i reaches di, and we set K{ = di. The 
trajectory of i is the segment sp^; in other words it is the track of i in the motorcycle graph. 

In the original motorcycle graph problem, the destination point di is at infinity in direction 
r>i. We can handle this case by computing a bounding box that includes all the vertices of the 
arrangement of the supporting lines £i, i = l,...,n, and choosing as destination points the 
intersections of the rays (sj, Vi) with the bounding box. The bounding box can be computed in 
0(n log n) time as any extreme vertex in the arrangement is the intersection of two lines with 
consecutive slopes. 

Unless specified otherwise, we make the following general position assumptions. No two mo- 
torcycles share the same supporting line, or have parallel supporting lines. No three supporting 
lines are concurrent. No point Sj, di lies on lj if j ^ i. No two motorcycles reach the same 
point at the same time. (We make these assumptions so as to simplify the description of the 
algorithm and the proofs, but our results still hold in degenerate cases.) 

The crossing point Xij is the intersection between ti and £j , and thus Xij = Xji =4n lj . 
The size \pq\ of a segment pq is the number of crossing points Xij that lie in pq. (See Figure [5j) 
We will need a data structure to answer halving queries: Given a query (i,p,q) where p, q are 
points on the supporting line £i, find a point h = h(p,q) 6 pq such that \ph\ ^ [pIjP?!! ana - 
\hq\ ^ [pIp^H, for a constant p < 1. In addition, we require that h is not a crossing point, and 
that both \ph\ and \hq\ are strictly smaller than \pq\ if \pq\ 2. 

2 Algorithm for Computing Motorcycle Graphs 

In this section, we present our algorithm for computing motorcycle graphs, as well as its proof of 
correctness and analysis. An example of the execution of this algorithm on a set of 4 motorcycles 
is given in Appendix [B| 

2.1 Algorithm description 

Our algorithm maintains, for each motorcycle i, a confirmed track SjCi, and a tentative track 
s{ti, such that Cj £ Sidi and tj G Cjdj. So the tentative track is at least as long as the confirmed 
track. As we will show in the next section, the confirmed track is a subset of the trajectory, 
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Figure 5: The size of pq is \pq\ = 3. Point h is a possible result h(p, q) of a halving query (5,p, q) 
with p = 1/2. 

so we have Ci 6 SjKj at any time during the execution of the algorithm. The tentative track, 



Our algorithm will ensure that no two tentative tracks cross. We keep all the tentative 
tracks in a ray shooting data structure, so that we can enforce this invariant by checking for 
intersection each time we try to extend a tentative track. This data structure returns the first 
tentative track hit by a query ray (p, v) , if any. We also build a data structure to answer halving 
queries, which will be used to shorten tentative tracks and keep them disjoint. 

Our algorithm builds the motorcycle graph by extending the confirmed tracks until they 
form the whole motorcycle graph. We may also update the tentative track of a motorcycle 
when we extend its confirmed track. 

A set of target points is associated with each motorcycle i. In particular, we maintain in 
a stack Si the set of target points that lie beyond the confirmed track of motorcycle i, thus 
Si C Cidi\{ci}. In other words, Si records the target points that motorcycle i has not reached 
yet. (See Figure[6j) The stack Si is ordered from a to dj. We denote by Top(Sj) its first element, 
so Top(5j) is the target point in Si that is closest to Cj. At the beginning, we set Si = {s{,di} 
for all i. New target points will be created in Case (3b) of our algorithm, as described below. 

If motorcycle i has neither crashed nor stopped, then its tentative track ends at the first 
target point in Si, so tj = Top(S'j). Otherwise, the tentative track and the confirmed track are 
the same, thus U = ty. So after a motorcycle has crashed or stopped, the ray shooting data 
structure records its confirmed track. 

An event (i,p) happens when a motorcycle % reaches a target point p. We process events one 
by one, and while an event is being processed, new events may be generated. After an event 
has been processed, we process the earliest available event. As t% = Top (Si) is the closest target 
point to % in Si, it means that we always process the event (i,U) such that r(i,ti) is smallest. 
Note that it does not imply that our simulation is done in chronological order: When we process 
an event we may create a new event (j,p) such that r(j,p) < r(i,ti). (See Appendix |Bj) 

We record in a priority queue Q the event (i, tj) for each motorcycle i that has not crashed 
or stopped. An event with earlier time r(i,ti) has higher priority. As ti = Top (Si), we can 
update the event queue Q in O(logra) time each time a stack Si is updated. So we can find the 
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Figure 6: This is the same example as Appendix|Bj Figure (m). The confirmed tracks are solid, 
and the tentative tracks are dashed. For motorcycle 1, the confirmed track and the tentative 
track go to X12 = ci = ti = Hi- The stack S± only records d\. For motorcycle 2, the confirmed 
track ends at C2 = Xi2> the tentative track ends at t2, and S2 = (^25X23,^2)- For motorcycle 
3, we have C3 = X23> £3 = ^3, and 53 = (d^). For motorcycle 4, we have C4 = S4, £4 = d^, and 
5*4 = (di). 

next available event in O(logn) time. The first n events are the events (i, Sj), i = 1, . . . , n, and 
occur at time t = 0. We process these n events in an arbitrary order. 

We now explain how to process an event (i, ti). To avoid confusion, for any motorcycle j, we 
use the notation Cj,tj to denote the endpoints of its confirmed and tentative track just before 
processing this event, and we use the notation c'j,t'j for their position just after processing this 
event. We first extend the confirmed track of motorcycle i to ti, thus = ti. We also delete U 
from Si. We are now in one of the following cases: 

(1) If ti = di, then motorcycle i stops. In order to avoid processing irrelevant events in the 
future, we remove Si from Q. 

(2) If ti is a crossing point Xij that lies in the confirmed track of j (that is, ti G SjCj), then i 
crashes at ti. So we remove Si from Q. 

(3) Otherwise, we try to extend the tentative track to the next target point q = Top(S'j). So 
we perform a ray shooting query with ray (ti,Vi), which gives us the first track intersected 
by Uq, if any. 

(3a) If tiq does not cross any track, then t[ = q, and we do not need to do anything else 
to handle this event. 

(3b) Otherwise, let j be the result of the ray-shooting query, so Sjtj is the first track hit 
by segment tiq, starting from ti. We shorten the tentative track of i, which means 
that we insert the new target point Xij into Si, as well as the point ^ = h(ti,Xij) 
obtained by a halving query on Uxij- If the crossing point Xij does not lie in the 
confirmed track of j, that is, if Xij £ c jtj \ i c j}i then we also shorten the tentative 
track of j, so we insert Xij into Sj, and we insert t'j = h(cj, Xij) into Sj. 

After applying the rules above, we update the ray shooting data structure (if needed), and we 
move to the next available event. 
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2.2 Proof of correctness 



Initially, we create the target points Sj, di for i = l,...,n. After this, we create new target 
points only in Case (3b) of our algorithm. There are two types of such target points: the 
crossing points Xij obtained by ray-shooting, and the points obtained by halving queries. We 
call x-targets the first type of target points, and /i-targets the latter. By our assumption that 
the result of a halving query is not a crossing point, a target point cannot be both a /i-target 
and a x-target. 

We need the following lemma. Remember that we say that two segments cross if their 
relative interiors intersect. 

Lemma 1. During the course of the algorithm, no two tentative tracks cross. 

Proof. For sake of contradiction, assume that two tentative tracks cross during the course of the 
algorithm. Let (i, ij) be the first event that generates such a crossing, so just before processing 
this event, the tentative tracks sdj, j = 1, . . . ,n do not cross, and there is a crossing among 
the tracks Sjt 1 -. We must be in Case (3), because we do not extend any tentative track in cases 
(1) and (2). Besides, we only extend the track of motorcycle i in Case (3). So there must be 
another motorcycle k ^ i such that Sji^ crosses s^tk- 

In Case (3a), the segment Uq obtained by ray shooting does not cross any tentative track 
Sjtj, j / i, and since t! i = q, then the new portion t^ of the track does not cross any other 
tentative track. The same is true in Case (3b), because t\ is in Uxij, where track j is the first 
track hit by tiq. So we just proved that, in any case, the new portion tit\ of the track does not 
cross any track Sjtj, j ^ i, and in particular, trf^ does not cross s^tk- 

By our assumption, we also know that Sktk cannot cross SiU. So the only remaining possi- 
bility is that s^k crosses at U. Then ti is the crossing point Xik- This point U = Xik cannot 
be in the confirmed track s^Cfe, because that would be Case (2) of our algorithm, and we showed 
that we are in Case (3). Since Xik is a x-t ar g e t of i, and it does not lie in the confirmed track 
of k, then it must have been inserted at the same time in Si and Sk while processing a previous 
event. Since Xik is not on the confirmed track of k, then it must still be in Sk- So the tentative 
track s^k cannot contain Xik i n hs relative interior, a contradiction. □ 

We want to argue that our algorithm computes the motorcycle graph correctly. So assume 
it is not the case. As our algorithm moves motorcycles forward until they either reach their 
destination point or crash, it could only fail if during the execution of our algorithm, the 
confirmed track of at least one motorcycle i goes beyond the point m where it is supposed to 
crash in the motorcycle graph. Let us consider the event (i,U) that is first processed by our 
algorithm, such that motorcycle i goes beyond K{. So Ki is in the segment Cjij \ {ti}. Let j 
denote the motorcycle that i crashes into, in the (correct) motorcycle graph, so Ki = Xij- 

When we process (i,ti), in the current graph constructed by our algorithm, motorcycle j 
cannot have reached dj, because it would mean that the tentative tracks Sjtj and Sjdj are 
crossing at Ki = Xij > which is impossible by Lemma [Tj 

We now rule out the case where, when our algorithm processes (i,ti), motorcycle j has al- 
ready crashed into some motorcycle k in the graph constructed by the algorithm. (See Figure[7|) 
For sake of contradiction, assume it did happen. 

• If we had i = k, then Xij would have been created as a x-target for j earlier. At this point, 
i had not gone past Xiji because (i,ti) is the first such event. As r(j,Xij) < i~(i,Xij), the 
algorithm would have moved j to Xij before i moves further, and thus j would not crash 
at Xij-, a contradiction. 



S 




Figure 7: The motorcycle graph (left) and an incorrect computation (right). 




™i Xij 



Figure 8: Proof of correctness, remaining case. 

• Thus we must have i 7^ k. As ti is beyond Ki = Xiji anci tentative tracks cannot cross, we 
must have Cj S SjXij- So j crashed into k at Xjk £ 'sJXij- As in the correct motorcycle 
graph, j does not crash into k, it means that the algorithm has already moved k past its 
(correct) crashing point, which contradicts our assumption that (i,ti) was the first such 
event. 

We just proved that j has not stopped or crashed when the algorithm processes event (i,ti), 
so at this point there should be an event (J,tj) in the queue. By Lemma [lj the tracks s{ti and 
Sjtj cannot cross, so we must have tj G cJkh- (See Figure|8j) It implies that r(j,tj) ^ r(j,Kj). 
But since i crashes into j in the (correct) motorcycle graph, we must have r(j, Kj) < r(i,Ki), 
thus r(j,tj) < r(i,Kj). As K{ G Cjij, we have r{i,Hi) ^ r(i, U), thus r(j,tj) < r(i,ti). But this 
is impossible, because our algorithm always processes the earliest available event, so it would 
have processed (j,tj) rather than (i,ti). 

2.3 Analysis 

Our algorithm uses two auxiliary data structures: for answering halving queries, and for ray 
shooting. The running time of our algorithm depends on their preprocessing time and query 
time. Let P(n) denote an upper bound on the preprocessing time of these two data structures, 
and let Q(n) denote an upper bound on the time needed for a query or update — so we can answer 
a ray-shooting query or a halving query in time Q(n), and we can update the ray shooting data 
structure in time Q(n). We now prove the following result: 

Theorem 2. We can compute a motorcycle graph of size n in time 0(P(n) + n(Q(n) + 
log n) log n) . 

Each time we handle an event, we perform at most two halving queries, one ray-shooting 
query, and we may update two tentative tracks in the ray-shooting data structure. We also pay 
an O(logn) time overhead to update the priority queue Q. So after preprocessing, the running 
time will be at most the number of events times Q(n) + logn. Thus we only need to argue 
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Figure 9: Proof of Lemma 3l first case. 



that our algorithm processes a total of 0{n log n) events. In fact, at each event we process, a 
motorcycle reaches a target point, so we only need to show that 0(n log n) target points are 
created during the course of the algorithm. 

Initially, we create 0(n) target points, which are Sj, di for i = 1, . . . , n. After this, we only 
create new target points in Case (3b) of the algorithm. In this case, we create one x-target, 
and at most two /i-targets obtained by halving. Thus we only need to bound the number of %- 
targets. At the end of the algorithm, some of these x-targets Xij correspond to an actual crash, 
with motorcycle i crashing into j, or j crashing into i. In any case, there are at most n such 
X-targets. We need to consider the other x-targets, that do not correspond to an actual crash. 
In this case, either motorcycle i or j does not reach Xij> so at the end of the computation, Xij 
must appear in the stack S{ or Sj of target points that have not been reached by motorcycle % or 
j, respectively. Thus, in order to complete the proof of Theorem [2j we only need the following 
lemma. 

Lemma 3. At the end of the execution of our algorithm, for any motorcycle i, the number of 
X-targets in Si is O(logn). 

Proof. In this proof, we only consider the status of the stack Si at the end of the algorithm, 
and we assume that it contains more than one x-target. "We denote by Xl, • • • , Xm the x-targets 
in Si, in reverse order, so Xm ■ ■ ■ X2X1 is a subsequence of Si, where Xm is closest to Ki and xi 
is closest to di. 

Each target Xj was created in case (3b) of our algorithm. At the same time, an /i-target 
hj = h(gj,Xj) was created by a halving query using another target point gj. As the points 
Xj, j = 1, . . . ,m are in Si, motorcycle i never reaches these points during the course of the 
algorithm, so xi an d hi must have been created first, then X2 an d hi . . . and finally Xm an d h m . 

For any 2 ^ j ^ m, as Xj 1S created after Xj-ii an d these two points are created when 
motorcycle i reaches gj and Qj-\, respectively, it implies that gj-i is in SigJ. We also know that 
Xj-i h es in Xjdi, because Xj-i appears after Xj i n Si. So gjxj, j = I, ■ ■ ■ ,m is & sequence of 
nested segments, that is, we have gjxj C gj-iXj-i for all 2 ^ j ^ m. More precisely: 



If hj-i is in Si, then gjXj C gj-\hj-i, because Xj is created after hj-\, and motorcycle i 
never reaches hj-\. (See Figure [9|) 



If hj-\ is not in Si, then gjXj C hj-ixj-i- (See Figure 10 ) It can be proved as follows 



As hj-% is created at the same time as Xj— 1> then Xj is created after hj-\. So Xj must 
have been created after motorcycle i reaches hj-t, otherwise we would have Xj £ Sihj-i, 
and since motorcycle i reaches hj-% later, Xj would not be in Si. As Xj is created after 
motorcycle i reaches hj—i, we must have gj £ hj-iXj-i- 



Thus gjXj is contained in either gj-\hj-\ or hj-iXj-lt an d since hj-\ = h{gj-\,Xj-i)i it 
follows that the size |<7jXjl decreases exponentially when j increases from 1 to m. As gm-iXm-i 
contains Xm a nd Xm-i, we have \g m -iXm-i\ j3 2. In addition, | <7iXi I ^ n i so we must have 
m = 0(log n). 

□ 
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Figure 10: Proof of Lemma 3l second case. 



3 Auxiliary Data Structures 



Our algorithm, presented in Section 2.1, requires two auxiliary structures. The first one is 
simply a ray-shooting data structure. As ray shooting is a standard operation in computational 
geometry, we will be able to directly use known data structures. The second data structure we 
need is for answering halving queries. We show below how to construct efficient data structures 
for this type of queries, and the corresponding time bounds for our motorcycle graph algorithm. 



3.1 General case 

In this section, we present the auxiliary data structures for the most general case, as presented 



in Section 1.2. So motorcycles have arbitrary starting position, destination point and velocity. 

For ray shooting, we can directly use a data structure by Agarwal and Matousek |2j, which 
requires preprocessing time 0(n 4//3+e ), with update and query time 0(n 1//3+e ), for any e > 0. 

For halving queries, we use known range searching data structures and parametric search, as 
in the work of Agarwal and Matousek on ray shooting: Our problem is an optimization version 
of range counting in an arrangement of lines, so we obtain the same bounds [2] Section 3.1]. 

Lemma 4. Given the n supporting lines i\,...,l n , we can construct a data structure with 
0(n 4/,3+e ) preprocessing time and 0{n l / i+£ ) query time that answers the following queries 
(i,p,q). Assume there are k crossing points Xij on PQ- Then we return the median cross- 
ing point and the next: the \k/2]th and the (\k/2~\ + \)th such crossing point, in the ordering 
from p to q along pq. 

With the two auxiliary data structures above, Theorem [2] yields the following result. 

Theorem 5. A motorcycle graph can be computed in 0(n 4//3+£ ) time, for any e > 0. 

It should be possible to replace the n £ factor in the bounds of Lemma|4]with a polylogarithm 
using known range searching techniques |11| [2"5] , because we only need a static data structure 
for halving queries, but in any case we need a dynamic data structure for ray shooting queries, 
so it would not improve our overall time bounds. 



3.2 C-Oriented Motorcycle Graphs 

We consider the special case where motorcycles can only take C different directions d\, . . . ,dc- 
Eppstein and Erickson gave an 0(ra 4 / 3+£ )-time algorithm when C = O(l). We show that 
with appropriate auxiliary data structures, we can solve this case in time 0(n log 3 n). In the 
following, we do not assume that C = 0(1), so our time bounds will also have a dependency on 
C. 

Proposition 6. We can compute a C -oriented motorcycle graph in 0(Cnlog 2 (n) min(C, logn)) 
time. 

We use the following data structures, and then the result follows from Theorem [2] 
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Ray shooting data structures. A first approach to answer our ray shooting queries is to use 
C instances of a data structure for vertical ray shooting in a planar subdivision. Several data 
structures are known for this problem [6], we use a data structure by Cheng and Janardan p3] 
that takes 0(log 2 n) time per update and O(logre) time per query. So overall, we get Q(n) = 
0(C log 2 n) with the terminology of Theorem^ 

Alternatively, we can use C(C — 1) instances of a data structure for vertical ray shooting 
among horizontal segments. Each data structure is used to answer ray shooting queries with a 
given direction, into segments with another direction: We just need to change the two coordinate 
axis to these two directions. Using a recent result by Giyora and Kaplan [22] . we obtain 
Q(n) = 0{C 2 log n). 



Halving queries. Our data structure for halving queries simply consists of a sorted list 
of motorcycles for each direction. So for each k 6 1,...,C, we have an array Ak of the 
motorcycles with direction dk, sorted according to the intercept of their supporting lines with 
a line orthogonal to dk- We now explain how to answer a halving query pq C li- 

Without loss of generality, assume ii has direction d\. For each direction dk, k = 2, . . . , C, 
the subset of motorcycles whose supporting lines cross pq appear in consecutive positions in 
Ak- We can find the first and the last index of these lines in O(logra) time by binary search. 
So we obtain all the arrangement vertices in pq in C — 1 sorted subarrays. We then compute 
the median rrtk of each such subarray Ak H pq, and the median of these points m& weighted by 
the number of points \Ak <^pq\ i n the corresponding subarray. This gives a halving point h(p, q) 
with p = 3/4. The median of each subarray can be found in O(l) time, and their weighted 
median in 0{C) time |17j . so the query time is dominated by the C binary searches. Thus, we 
can answer halving queries in 0{C log n) time. 



3.3 Bounded precision input 



The data structure for answering halving queries in Section 3.1 is quite involved. In practice, 
one would rather implement halving queries by simply halving the Euclidean length \\pq\\ instead 
of approximately halving the number of crossing points. Unfortunately, in the infinite precision 
model that is commonly used in computational geometry, this would cause the analysis of our 
algorithm in Lemma [3] to break down, because a stack of target points Si may have size Q{n) 
at the end of the algorithm. 

Such a counterexample would require the distance between consecutive target points in Si 
to become exponentially small near the crashing point, which does not seem likely to happen in 
practice. To formalize this idea, we make the assumption that all input numbers (the coordinates 
of the starting points, the destination points, and the velocities) are rational numbers, whose 
numerator and denominator are in {— 2 W ~ 1 , . . . , 2 W ~ 1 — 1} for some integer w. In other words, 
the input numbers are w-bit signed integers. We still assume that arithmetic operations between 
two numbers can be performed in constant time. 

This model also allows us to handle the case where input coordinates are u>-bit rational 
numbers, that is, rational numbers with w-h\t numerator and denominator; we just need to 
scale up each coordinate by a factor 2 W to obtain 2w-bit integers, losing only a constant factor 
in our time bounds. In the proofs below we assume the input numbers are integers, to simplify 
the presentation, but the results are stated with rational coordinates. 

As the input coordinates are w-bit integers, the coordinates of a crossing point Xij are ratio- 
nal numbers obtained by solving a 2 x 2 linear system, their denominator being the determinant 
det(vi, Vj). Thus, the denominator is an integer between — 2 2w ~ 1 and 2 2w ~ 1 . So any two distinct 
crossing points are at distance at least 2~ 2w+l from each other. 
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Assume that we replace our halving operation, as defined in Section 1.2 with halving the Eu- 
clidean length. So h{p, q) is the midpoint oipq, which can be computed in constant time. Then 
any nested sequence of segments obtained by successive halving, as in the proof of Lemma [3j 
consists of 0(w) such nested segments, because a segment g~Jx] °f length smaller than 2~ 2w+1 
cannot contain another crossing point in its interior, and hence it cannot be further subdivided. 
So the bound on the size of Si becomes 0(w), and we get the following result. 

Theorem 7. If the input coordinates are w-bit rational numbers, we can compute a motorcycle 
graph in time 0(nw(Q' (n) +logn)), where Q'(n) is the time needed for updates or queries in 
the ray-shooting data structure. 

For bounded-precision input, the bottleneck of our algorithm has thus become the ray shoot- 
ing data structure, whose update and query time bound is 0{n l ' i+£ ) in the most general case. 
Therefore, we obtain a faster motorcycle graph algorithm if we are in a special case where faster 
ray-shooting data structures are known. One such case is ray-shooting in a connected planar 
subdivision, which can be done in 0(log 2 n)-time per update and query using a data structure 
by Goodrich and Tamassia [23] . We can use this data structure if, for instance, all motorcycles 
move inside a simple polygon P, starting from its boundary. (So for all i, sidi C P, and Sj is 
on the boundary of P.) Then we perform ray shooting in the union of the tentative tracks and 
the edges of P, which form a connected subdivision. It yields the following time bound. 

Corollary 8. We can compute a motorcycle graph in time 0(nlog 3 n) for n motorcycles mov- 
ing inside a simple polygon with 0{n) vertices, starting on its boundary, and if the input has 
0(\ogn)-bit rational coordinates. 



4 Application to Straight Skeleton Computation 

In this section, we give new results on straight skeleton computation, using our new motorcycle 
graph algorithm. 



4.1 Preliminaries and non-degenerate cases 

As we mentioned in the introduction, the straight skeleton problem and the motorcycle graph 
problem are closely related. We now explain it in more details. 

Consider the reflex (non-convex) vertices of a polygon P. When we construct the straight 
skeleton of P, these vertices move inward and may collide into edges, or other vertices. These 
events, called split events and vertex events, are the difficult part of straight skeleton compu- 
tation, because they affect the topology of the shrinking polygon by splitting it, and because 
they are non-local: A reflex vertex may affect a chain of edges on the other side of the poly- 
gon. The other type of events, called edge events, where an edge shrinks to a point, are easily 
handled with a priority queue. So the interaction between reflex vertices is a crucial part in 
straight skeleton computation, and the motorcycle graph presented below helps to determine 
these interactions. 

The motorcycle graph induced by a polygon P is such that each motorcycle starts at a reflex 
vertex, moves as the same velocity as the corresponding reflex vertex when we shrink P, and 
stops if it reaches the boundary dP of P. (See Figure [2j) 

If P is degenerate, then two reflex vertices may collide and create a new reflex vertex. In 
this case we need to create a new motorcycle after the collision. (See Figure |3j) So when two 
motorcycles collide in the induced motorcycle graph, we may have to create a new motorcy- 
cle [25]. Our motorcycle graph algorithm, as described above, does not apply directly to this 
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Figure 11: An example where S% holds a linear number of target points at the end of execution 
of the algorithm. The speed of the motorcycle at the bottom are adjusted so that xi is created 
first, then X2 - ■ ■ 



case, because the proof of Lemma [3] breaks down. (The reason is that Si may hold a linear 
number of target points at the end of the execution of the algorithm, due to the newly created 



motorcycles. See Figure 11.) In Section 4.2, we will explain how to compute these generalized 
motorcycle graphs efficiently on bounded-precision input. But the following theorem still holds 
in degenerate cases. 

Theorem 9 (Cheng and Vigneron [15j . Huber and Held [25]). The straight skeleton of a polygon 
P with n vertices and h holes can be computed in expected time 0(n\Jh + 1 log 2 n) if we know 
the motorcycle graph induced by the vertices of P. 

From the discussion above, and using our motorcycle graph algorithm from Theorem [5j we 
obtain the following result. 

Corollary 10. We can compute the straight skeleton of a non-degenerate polygon with n vertices 
and h holes in 0(n 4//3+£ + n\/h + 1 log 2 n) time for any e > 0. 

4.2 Bounded precision input 



We use the same bounded precision assumptions as in Section 3^3 where the input coordi- 
nates are w-bit integers or, equivalently, ui-bit rational numbers. Similarly, to simplify the 
presentation, we use the integer model in the proofs, but we state the results in the rational 
model. 

Thus, the coordinates of the vertices of the input polygon P are w-bit integers. The sup- 
porting lines li of the motorcycles are angle bisectors between two edges of the input polygon. 



In order to apply the same halving scheme as in Section 3.3, where the Euclidean length is used 
instead of the number of arrangement vertices, we need to argue that the separation between 
two vertices in this arrangement of bisectors cannot be too small. This distance can be shown 
to be at least 2~ w , where W = 64(80u> + 105) + 1 = 0(w), by applying the separation bound 
by Burnikel et al. [9|. So we obtain a result for induced motorcycle graphs that is analogous to 
Theorem 

Lemma 11. Given a polygon P whose input coordinates are w-bit rational numbers, we can 
compute the motorcycle graph induced by P in time 0(nw(Q'(n) + logn)), where Q'(n) is the 
time needed for updates or queries in the ray-shooting data structure. 

In the lemma statement above, we do not exclude degenerate cases. This is another ad- 
vantage of this bounded precision model. As the argument in our analysis only relies on the 
separation bound between two distinct crossing points, and not on the number of motorcycles 
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crossing a given segment, a newly created motorcycle does not affect our analysis as it still 
obeys the same separation bound: A newly created motorcycle still lies on the bisector of two 
input edges, though these two edges are not adjacent in the input polygon [25]. (See Figure [3}) 

We still need to describe an efficient ray-shooting data structure. As our input polygon has 
h holes, the boundary dP of P together with the tentative tracks form a collection of h + 1 
disjoint simple polygons. We could directly use known ray-shooting data structures [U |2"4"] . 
which can be made dynamic at the expense of an extra n e factor in the running time pQ. In 
the following, we give a different approach, which leads to a better time bound when used as 
a subroutine of our algorithm. This approach takes advantage of the fact that the holes of P 
are fixed (only the tentative tracks are dynamic). We use a spanning tree with low crossing 
number, which is not a new idea in ray-shooting data structures |13| 124] . 

We pick one point on the boundary of each hole of P, and on the boundary of P. We 
connect these h + 1 points using a spanning tree T with low stabbing number [3J, that is, a 
spanning tree such that any line crosses at most 0{s/h) edges of T. This tree can be computed 
in 0{n 1+£ ) time [3J Section 8]. We maintain a polygonal subdivision which is the overlay of P 
with T and the tentative tracks. So at each intersection between an edge of T and an edge of 
P or a tentative track, we split the corresponding edges and tracks at the intersection point. 
This subdivision S is connected and has 0(nVh) edges, and we maintain it in the ray shooting 
data structure by Goodrich and Tamassia [23], which has 0(log 2 n) update and query time. 

Each time a tentative track is extended or retracted, as a tentative track intersects O(Vh) 
edge of T, we can update the subdivision and the data structure by making 0{\/h) updates in 
the ray shooting data structure. Similarly, when our motorcycle graph algorithm tries to extend 
a tentative track, we can find the first tentative track being hit by a query ray in 0(Vhlog 2 n) 
time: We first perform a ray shooting query in S, which takes 0(log 2 n) time. If we hit an edge 
of T, we make another ray shooting query starting at the hitting point of the previous query, 
and in the same direction. We repeat this process as long as the result of the query is an edge 
of T , and by the low-stabbing number property, it may only happen 0(y/h) times. 

Overall, our ray shooting data structure has update and query time 0{\fh log 2 n). So by 
Theorem [9] and Lemma 11 , we obtain the following result. Note that it still holds for degenerate 
input. 

Theorem 12. The straight skeleton of a polygon with n vertices and h holes, whose coordinates 
are 0(logn)-bit rational numbers, can be computed in 0(n\/h + 1 log 3 n) expected time. 
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Appendix 



A Pseudocode 

In this section, we give the pseudocode of our algorithm. It is more detailed than the algorithm 
description in Section [2. 1[ and it can handle degenerate cases. The proof of correctness and the 
analysis are essentially the same as in Section [2j but they require a more detailed case analysis. 

To deal with the degenerate cases where some supporting lines are concurrent, or two or 
more motorcycles reach a point at the same time, we record all the target points created so far 
in a dictionary data structure T>. We can implement T> as a balanced binary search tree, sorted 
in lexicographical order of the coordinates (x,y), which allows to retrieve a point in O(logn) 
time. We associate three fields with each point p stored in T>: 



A list M(p) records the motorcycles i such that p G Si. So M(p) records all the motorcycles 
i that could possibly reach p, at a given point of the execution of our algorithm. The set 
M(p) itself is stored in a dictionary data structure, so that we can decide in O(logn) time 
whether a motorcycle i is in M(p). 

• Two motorcycles k,k' of M(p) such that T(k,p),T(k' ,p) are smallest. It will allow us to 
find out whether two motorcycles crash simultaneously at p. 

• A flag Blocked (p) which is set to FALSE initially, and is set to TRUE as soon as a 
confirmed track has reached p, implying that any other motorcycle that reaches p must 
crash. 

After the initialization stage, our algorithm handles repeatedly the earliest available event, 



according to the four cases (1), (2), (3a) and (3b) from Section 2.1 



Lines 11 and 12 deal with Case (1) and (2). The condition p = d. L corresponds to Case 
(1). The other two conditions check whether we are in Case (2). In particular, condition 
Blocked(p) = TRUE means that at least one other motorcycle has reached p, thus motorcycle i 
crashes. With degenerate input, it is possible that another (or several other) motorcycle reaches 
p at the same time as i, in which case Blocked(p) = FALSE if (i,p) is the first event involving 
p that has been processed. The condition at Line |12| checks whether we are in this situation. If 
so, i must crash. As (i,p) is the first event involving p that we process, there is no earlier event 
in M(p), so we can find another motorcycle that reaches p at the same time in constant time 
using the second field associated with p in V. 



Case (3a) corresponds to a positive answer to the test at Line 19. The condition d(si,p') > 
d(si, Top(Si)) detects whether the track of i to the next target points hits any other track. The 
other condition p' = Top(S'j) and j G M(p'), checks for a boundary case, where Top(5j) falls on 
another tentative track. The test is positive if p' has already been identified as a target point 
of j before. In this case we only extend the tentative track of i, without doing any unnecessary 
halving. 



Line 21 branches to Case (3b). Similar to Line 19, we do not perform an unnecessary halving 
operation when i G M(p') or j G M(p'). 

Our pseudocode does not handle explicitly the case where two motorcycles have same sup- 
porting line. These cases can be easily handled by ad-hoc arguments |15j . One way of doing it 
is to insert additional target points at initialization. For each supporting line shared by several 
motorcycles, between any two consecutive motorcycles i,j along this line that go toward each 
other, we insert their potential collision point, that is, we insert into Si and Sj the point p such 
that r(i,p) = r(j,p). For each motorcycle i along this line, if the first starting point Sj in the 
ray (si,Vi) is in Sidi, we also update di to be Sj. 
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Algorithm 1 motorcycle_graph 

Initialize the dictionary D. > Initialization 

for i = 1 — > n do 

Set q <— Si, ti <— Si and Si <— {si,di}. 
Insert Si and di into X>; 
Insert motorcycle i into M(sj) and M{di). 

Initialize the event queue Q and the data structures for ray-shooting and halving, 
while Q is not empty do > Main loop 

Let (i,p) be the earliest available event. > So p = ti. 

Set Cj <— p. 

Pop Top(Si) from > Here Top(^) = p. 

if Blocked(p) = TRUE, or p = dj, or 

3k G M(p) \ {z} such that r(k,p) = r(i,p) then 
Motorcycle i crashes. 
Set ti ^— p. 

Remove i from M(j) for all j G 5j. 
Remove Sj from Q. 
else 

<- rayshooting(cj,-yi). 
if d(si,j/) > d(ai,Top(Si)), or (p' = Top(Si) and j G M(p')) then 

Set ^ <- Top(5j). 
else 

if i G M(p') then 

Set ^ <- Top(5i). 
else 

Push p' into Si . 
Push p* = h(ci,p') into 
Set ti ^— p* . 

Insert i into M(p') and M(p*). 

if p' ^ s^cj and j ^ M(p') then 
Push p' into Sj. 
Push p* = h(cj,p') into 5,-. 
Set tj <— p* . 

Insert j into M(p') and M(p*). 
Set Blocked(p) <- TRUE. 
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B Example 



We give an example of the execution of our algorithm on a set of 4 motorcycles. (Confirmed 
tracks are solid, and tentative tracks are dotted.) It demonstrates two features of our algorithm, 
that were mentioned above. 

• A tentative track may be longer than the final track in the motorcycle graph. For instance, 
the tentative track s\di in (b) is longer than the final track S1X12 in (q). 

• Our algorithm does not construct the motorcycle graph in chronological order. For in- 
stance, in (i), motorcycle 2 is moved to X12, which is its position at time r(2,xi2) = 
2.12072. Then in (k), motorcycle 3 is moved to p4, which is its position at time t(3,£>4) = 
1.667206. 

The four motorcycles 1,2,3 and 4 start at time at initial points si = (0.8,3.3), S2 = 
(0.5,1), S3 = (5.7,0) and S4 = (6,3.4). Their velocities are v± = 1.2(cos — 5°, sin — 5°), V2 = 
1.7(cos35°,sin35°), v 3 = 2(cos93°,sin93°), and v 4 = 0.8(cos -37°, sin -37°). 

We use the halving scheme as specified in Section 1.2 with p = 1/2. So for instance, we 
create p 4 in (j) by halving 53X23- There are three crossings along this segment: Xi3>X23,X34- 
Then p 4 is created as a point between X13 and X34> m this case we just use the midpoint. 
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r(2,xi 2 ) = 2.120721 



(d) new events: r(l,pi) 
r(2,p 2 ) = 1.060361 



= 1.109735, and 
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dz d 2 



d-t d 2 



Sl 


Pi 


Xl2 
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S*4 
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•■' 








S2 









S3 

(e) r(3,d 3 ) = 2.503431 
<^3 d 2 



(g) move: r(2,p 2 ) = 1.060361 

<^3 <^2 
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Pi 


XV2 \ 
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S*4"... 
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• X23 


•Sl 


Pi 


Xi?. ■■- 


; s'i-.... 




yV2 






S2 









C?4 



(i) move: r(2,xi2) — 2.120721. new events: 
t(2,x 23 ) = 3.565653, and t(3,x 23 ) = 2.24 147 



• P. 1 X12 

Si 


S*4" 


• 

P2 




* 

S2 






53 


(f) r(4,d 4 ) = 3.130339 

d 3 C?2 


. Pi X12 

Sl - 


• 

si--,. 


^^^/¥2 




S2 




S3 


(h) move: r(l,px) = 1.109735 
d 3 d2 


• 


X2 3 


• P. 1 Xl2 Pi • 


S*4-,. 










S3 


(j) new events: r(2,p 3 ) = 


2.676739, a 



r(3,p 4 ) = 1.667206 
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d 3 d 2 rf 3 d 2 




S3 

(q) move: r(2,x 2 3) = 3.56 5 6 53 
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